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Abstract 

We compute the Schrodinger functional (SF) for the case of 
lattice QCD with Wilson fermions (with and without SW im- 
provement) at two-loop order in lattice perturbation theory. This 
allows us to extract the three-loop /3-function in the SF-scheme. 
These results are required to compute the running coupling, the 
A-parameter and quark masses by finite size techniques with neg- 
ligible systematic errors. In addition our results enable the im- 
plementation of two- loop 0(a) improvement in SF-simulations. 
This article is based on the revised version of ref . [1 1] . 



1 Introduction 



Physical amplitudes in quantum cliromodynamics (QCD) obey renormal- 
ization group equations whose solutions at high energy can be computed (in 
various leading logarithmic approximations) using renormalized perturbation 
theory. The solutions involve renormalization group invariant mass parame- 
ters, the so-called (scheme dependent) A parameter and the renormalization 
group invariant quark masses M. The determination of the relationships of 
these quantities to the low-lying hadron masses is a non-perturbative prob- 
lem. For the quark masses there are various non-perturbative approaches 
including chiral perturbation theory, but the most systematic approach and 
the only one available for computing e.g. the ratio A/Mproton is to employ 
the lattice regularization. There is already a vast literature on this subject 
(see e.g. the reviews [fl-p)- 

Among the collaborations which have as one of their main goals the com- 
putation of these relationships to as high precision as presently possible is the 
Alpha Collaboration. The collaboration's results in the quenched approxi- 
mation have already been presented, most recently in [Q. The extension of 
the project to full QCD is now under way but due to the enormous extra 
cost in CPU time, without the sufficient compensation of computer power it 
will still take some time until one can match the accuracy attained in the 
quenched approximation. 

The emphasis of the Alpha project is on precision and directly associ- 
ated with this, an attempt to carefully control systematic errors. The aim is 
to compute non-perturbatively defined running couphngs and quark masses 
over a wide range of energies ranging from low energies where contact with 
hadronic mass parameters is made. To make contact with perturbation the- 
ory it is essential that the simulation reaches a range of high energy where 
the predicted behavior actually appears to set in. Our tactic to reach this 
goal has been explained in many publications , [§] and we refer the reader 
to these papers for further technical details. The key idea is to define a cou- 
pling g{L) in a system with finite linear extension L (with specified boundary 
conditions). Adjusting a/L and the bare coupling so that g{L) is kept fixed, 
the approach to the continuum limit of the coupling at twice the extension 
g{2L) depends only on a/L. As a working hypothesis this approach is fitted 
by a power law (up to logarithmic factors) as expected from the Symanzik 
effective action analysis 

A specific such coupling osf = ^7sF/4vr, measured by the Alpha Collab- 
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oration is not directly experimentally accessible, but it can be related to 
phenomenological running couplings e.g. the MS coupling of dimensional 
regularization defined in infinite volume at high energies, using renormalized 
perturbation theory: 

^Wsi^^) = "SF(g) + ci(s)asF(g)^ + C2{s)asF{qf + ■■■ (1-1) 

To reduce the estimated systematic errors introduced by truncation of the 
perturbative series to the level of ~ 1% it is necessary to work out this 
connection to two-loop order. For technical reasons this is at present accom- 
plished in two steps. The first step, the computation of the relationship of 



a- 



MS 



to the bare lattice coupling ao, has been completed in refs. and 
ri|. In this paper we perform the second step, the relationship of asF to 



for A^f > flavors (the case A^f = has already been presented in |T^-[|I^). 

The final result between the continuum quantities (|l.lj ) is of course regu- 
larization independent, but the intermediate relations in both steps depend 
on the details of the lattice action employed. In our simulations we attempt to 
reduce lattice artifacts in the continuum limit by working with the Symanzik 
0(a)-improved fermion action of Sheikholeslami and Wohlert The can- 
cellation of extra 0(a) effects due to the presence of boundaries requires the 
introduction of additional boundary terms in the action. The weights of 
these terms can also be computed in perturbation theory and the computa- 
tion presented here fixes the one relevant for the boundary conditions under 
consideration to two-loop order. 

After recalling the definition of asp in the next section we will present 
some technical details, describing the Feynman diagrams involved in section 
3. In section 4 we discuss the dependence of the perturbative coefficients 
on the lattice cutoff, and the related task of determining the boundary im- 
provement coefficients mentioned above. In the final section 5 we compute 
the coefficients in eq. ( |1 . l\j . A summary of part of this work has already been 
given in 



16 . 



2 Basic definitions 

The Schrodinger functional is the QCD partition function 

Z = e-^ = J D[Um]D[ij]e-'' , (2.1) 
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with the following particular geometry. In the spatial directions we have a 
finite box of size L (in all directions) with periodic-like boundary conditions 
(see below). The extent of the time (xq) direction is also finite, T, and at the 
boundaries xq = 0, = T, Dirichlet boundary conditions (specified below) 
are imposed. The ratio T/L belongs to the definition of the SF and in this 
work we always restrict ourselves to the choice T = L. 

The SF can be straightforwardly defined non-perturbatively by using the 
lattice regularization. In eq. (|2.1|) the links interior to the box are integrated 
over with the invariant SU(3) measure, and the dynamical Grassmann fields 
ip{x),ilj{x) are those at points x with < xq < L. 

The action is taken to be 

S[U, ^, ^] = S[U] + Sw [U, ^, ^] + Ssw [U, ^, ^] . (2.2) 
The pure gauge part is defined by the usual sum over oriented plaquettes, 

SiU) = \Y.^{p)tT{l-U,). (2.3) 
9o p 

The weight w{p) is unity for all plaquettes except those at the boundary that 
contain the time-direction and one of the frozen spatial links where we put 

wip)=Ctigo). (2.4) 

S'w has the form of the standard Wilson fermion action with r = 1, 

^w[t/,^,^]= E ^{x){D + Mo{x))^{x), (2.5) 

x,0<xo<L-l 

where Mo(x) is the usual bare mass mo, modified only by extra contributions 
at the boundaries 

Mo{x) = mo + (^(5-0) - 1) [5x'o,i + Sxo,L-i] ■ (2.6) 
The lattice Wilson Dirac operator D is given by 

^ = ^E{7.(v; + v,)-v;v,}. (2.7) 

Hermitian 7-matrices are used (we adopt the conventions of ref. |l7]). V^, V* 
are the forward and backward lattice covariant derivatives 

V^^ix) =A^f/(x,/i)^(x + /i)-^(x), (2.8) 
V;ijix) =V>(x)-A*f/(x-/i,/i)V(x-/i). (2.9) 
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Note that to be able to write the action ( [2 .51 ) in a more elegant form, we 
have defined ip{x) = = ip{x) for xq = 0,L. The additional U{1) field in 

mm 

r 1 if u = 0, 

corresponds to a free phase in the spatial boundary conditions of the fermions. 
Finally Ssw is the Sheikholeslami-Wohlert |T^ term 



_ _ 



(2.11) 



where we refer to ref. ^Tj for the definition of the field J-'^^, appearing here. 
The function Csw(5'o) is assumed to have a perturbative expansion of the form 



.(0) 



+ cl^gl + ■■■ 



(2.12) 



Appropriate specification of Csw(5'o) allows the elimination of all 0(a) effects 
in on shell quantities; the values of the coefficients required for this improve- 
ment will be denoted by 42* • ^^^e level c^^* = 1 (or = r with general 
Wilson parameter). The value of the 1-loop coefficient 4^^* is known (see Ta- 
ble P and it is independent of Nf] we note also that c*^ has been computed 
non-perturbatively for a range of bare coupling qq in the quenched theory 
[]TB[, |T^ and in the unquenched theory with A^f = 2 |2C 



To define renormalized quark masses from the bare mass rriQ for Wilson 
fermions requires both additive and multiplicative renormalization: 



where the critical mass rric has a perturbative expansion: 



(2.13) 



rrir 



(2.14) 



The coefficients m^*-* depend on Cg^'i in particular the 1-loop coefficient m 



(1) 



.(0) 



'see Table 



which appears in our computation, depends on the tree value 4 
but is independent of Af. 

It remains to specify the boundary conditions. Firstly the spatial links 
U{x, k), k = 1,2, 3, are frozen to fixed SU(3)-elements for the layers at time 
coordinate Xq = 0,L, 



U{X, fc) 1x0=0 

U{x,k)\^^=L 



exp(mC) 
exp(zaC") 



exp 



exp 



la 



diag( 



^ diag( 



(2.15) 
(2.16) 
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coefficient 


value 


references 


'-SW ^ 


-0.4342856(3) 
-0.2700753495(2) 


24 


22 

1 1 


M 

25 


23| 

1 




0.26590(7) 


|24 


M 


25 


1 



Table 1: Perturbative coefficients for gauge group SU(3) 



Here a is the lattice spacing, and (pi, are certain p] dimensionless numbers 
that depend on one free parameter rj. The SF-coupling g is defined from the 
response in the free energy F ( |2.1| ), to infinitesimal changes in the surface 
fields by varying rj, 

f = k/V\ (2.17) 

where V is the derivative with respect to ?7 at 77 = 0, and is a constant 
which is fixed by normalizing the leading term in the perturbative expansion 



g\L)=g',+p,{I)gl + p,{I)gl + ... 



(2.18) 



where here and in the following we often denote the number of lattice points 
(in a given direction) L/a by /. 

The coefficients Cj(s) in the continuum relation (|1 . 1| ) depend not only 
on the number of quark fiavors A^f but also on the particular choice for 
the boundary phases 0,0' occurring in ( 2.15 . 2?T6| ) and on the parameter 9 
( |2.10| ). Both our perturbative calculation here and simulations reported in 
§ and more recently in P, |, |^ are restricted to the choice "A" of [H for 
the background field and 6 = n/5. 

The coefficients Pi appearing in (|2.18|) depend on the number of fiavors. 



Pi = Pw+PuNf, 

P2 = P20 + P2lNl + P22Nf'^ 



(2.19) 
(2.20) 



The one- and two-loop coefficients pio,P2o for the quenched case have been 
computed in refs. 0, [1^ respectively (and for the case of SU(2) in [p^] ); 
we note that in these papers the coefficients piQ were denoted by rrii. The 



one-loop coefficient pn was computed in ref. |]T^. Our objective here is to 
compute the remaining two-loop coefficients P2i,P22- 
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The Callan-Symanzik equation 



Lg^9{L) = -Pig) = bof + hf + hf + ... (2.21) 

(the universal 1- and 2-loop beta-function coefficients are given in Appendix 
A) requires for all values of the improvement coefficients that 

pi = 26oln(/)+pi + 0(l//), (2.22) 
P2-pI = 26iln(/)+p2 + 0(l//). (2.23) 



where pi are independent oi L/a but dependent on A^f and on Csw The 
known values including our results^ here are summarized in Table 0, where 
components pij are defined in analogy to ( |2.19| ), (|2.20|) . 



(^,J) 




<--sw — U 


P^3 


c<"'-l 


references 


(1,0) 

(1,1) 

(2,0) 

(2,1) 
(2,2) 


0.36828215(13) 
-0.009868186(4) 

0.048091(2) 
-0.00349(3) 

0.000211(3) 


0.36828215(13) 
-0.034664940(4) 
0.048091(2) 

0.00978(2) - 0.054641(l)cW 
0.000209(1) 


i 

|17] 
14 



Table 2: SF couphng coefficients for gauge group SU(3), background field 
A and = n/5 



The boundary weight Ct is assumed to have a perturbative expansion of 
the form 

c,{go) = l + ci'^g', + ci'^g', + ... (2.24) 

and similarly for Ct. The freedom of adjusting the boundary weights Ct and 
Ct to specific functions c* and c* is required for improvement of 0(a) lattice 
artifacts that are otherwise introduced by the surfaces. The one-loop coef- 
ficient c[^^* required for improvement of SF correlation functions has been 
computed in ref. |^ and is independent of A^f. With this knowledge, the 



perturbative coefficients of can be computed to 2-loops from evaluating 



^Reanalysis of the 2-loop Nf — data with the method of appendix D has led to errors 
smaller than the very conservative estimate of ||l4|. 
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the 0(a) lattice artifacts of the SF couphng and demanding deviations from 
leading continuum behavior to be of O(a^), 



c=c* 



26oln(/)+pi + 0(l//=^), 
26i ln(/)+p2 + 0(l//2). 



P2-P1 

The coefficients ^* do depend on the number of fermions, 



.(1)* 

-,(2)* 



^(1.0)* 



f ) 



(2,0)* I {2,1)* f,T I (2,2)*,,r2 



The values^] of these coefficients are given in Table ^. 



coefficient 


value 


references 


~(i)* 

Ct 


-0.01795(2) 


[25|,|27] 


H 
Ct 
H 

(2,1)* 
Ct 

J2,2)* 
Ct 


-0.08900(5) 
0.0191410(1) 

-0.0294(3) 
0.002(1) 
0.0000(1) 


i 



[14] 



(2.25) 
(2.26) 



(2.27) 
(2.28) 



Table 3: SF boundary 0(a) improvement coefficients for gauge group SU(3) 



3 Perturbation expansion to two loops 

The perturbative expansion in the bare coupling amounts to an expansion of 
the gauge field around the induced background field. This has been described 



in detail in previous publications [0], p6|, [T^. In particular we use the gauge 



fixing procedure as in ref. [0. In the presence of the abelian background field, 
both the gluon and quark propagators cannot be computed analytically. The 
numerical computation of the gluon propagator is described in Section 2 of 



^See previous footnote. 
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ref. [^1, and here in Appendix B we give some technical details on how we 
numerically computed the fermion propagator S. 

An enumeration of the Feynman diagrams contributing to the Schrodinger 
functional to two-loop order reveals that the coefficients pij depend on the 
improvement coefficients and critical mass in the following way: 



Pio 
Pu 

P20 - pIo 
P21 - 2pioPii 



Pi + ci'^'^pl , 



P2 



c, 



Mb 



P2 + 



ci 



(1,0) 



e , (1,1) b , o (1,0) (1,1) c , (2,1) d 



P22 



P?1 



P'2 + 



.(1,1) 



(3.1) 
(3.2) 

(3.3) 

(3.4) 
(3.5) 



All the coefficients coming from diagrams not involving the quarks have 
been computed in W^. Here we recall that only two of the coefficients can 
be given in closed form: 



P\ 



P2 



pi 



2 
7 



y-^ + 0(l//^) 



(3.6) 

(3.7) 



whereas the others require careful numerical evaluation. The contributions 
p1, P2 and P2 are tabulated in Table 1 of ref. [Q. 

As for the diagrams involving the quarks, the one-loop coefficient 



Pi 



^ -pel 

k ' 



with 



-pel 

^ 1 



-tr S 



dS- 
dr] 



1 77=0 , 



(3.8) 



(3.9) 



has been computed in ref. [17 



It remains to compute the coefficients pi, P2, P2j P2J P2^ pi- The contribu- 
tions P2,P2iP2 ^1'^ very simple because they are obtained by differentiating 
the one-loop contribution with respect to the bare parameters e.g. 



tr 



S——S 



S- 



dr] 9cs„(0) dr]dcs„{0) 



\ri=0 , 



(3.10) 
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and similarly for P2)P2- The contributions are diagrammatically represented 
in Fig. la where the filled circle indicates the insertion obtained by differen- 
tiation with respect to one of the bare parameters (the differentiation with 
respect to rj is not indicated). 

The contribution arises from the correlation of a boundary term in the 
pure gluon part of the action with the quark-antiquark 1-gluon vertex given 
in diagram Fig. lb. 




Fig. la Fig. lb 

The remaining part p| of p2i comes from the four diagrams in Fig. 2. The 
fermion big-mac diagram Fig. 2c is the technically most difficult diagram; 
the rest of the diagrams are simpler because they are essentially products of 
1-loop diagrams. 




Fig. 2c Fig. 2d 

Finally there is only one as yet unmentioned contribution P2 to P22, which 
is basically the product of two fermion tadpoles depicted in Fig. 3. 



9 



Fig. 3 



Despite optimization of the code in various respects by e.g. making use 
of all symmetry properties, the CPU-time for the computation which is dom- 
inated by the big-mac diagram in Fig. 2c, is rather costly growing rapidly 
as (L/a)^. All diagrams were computed using double precision arithmetic to 
at least L/a = 32, for two cases = 1 and = 0. The numerical results 
for pI,P2,P2,P2,P2,P2 for 4^ = 1, and for pl,pl,pi for 4^ = 0, are given in 
Appendix C in Tables and ^ respectively. 

We have various stringent tests (as in [^) to check our numbers; firstly 
the full code using sums in position space was compared on the smallest 
lattices with a slow but independent code (many of whose subroutines had 
been well tested in previous computations) making sums in momentum space. 
Secondly gauge parameter independence of the was checked. Thirdly it 
was checked that the full code reproduced numbers for partial sums satisfying 
all analytically known symmetries, before these were actually used to reduce 
the number of terms in order to speed up the program execution time. Finally 
the fact that the /-dependence for P2i,P22 (depending on the value of 4^ as 
described in the next section) was as expected from general considerations, 
gave a further consistency check. 



4 Analysis of the a/L dependence 

The generically expected behavior for 2-loop lattice Feynman diagrams sug- 
gests an asymptotic expansion for the two-loop coefficients pf of the form 
(recall / = L/a) 

CO 

P2 = E K + < ln(/) + t^n ln^(/)]/-". (4.1) 

The coefficients of these expansions are extracted from the series of finite 
lattices (see Appendix C). The method we used here to extract these numbers 
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together with an estimate of their systematic errors is described in some detail 
in Appendix D. 



4.1 The case 4^ = 1 

4.1.1 The continuum behavior of P22 

Let us start with the simpler case of P22 where the non-trivial contributions 
come from only and p{. For p^ we have ri'^ = 1 and t{ = for all n. Our 
result is: 

s{ = -0.03375(4), (4.2) 
r( = -0.0734(4). (4.3) 

The contribution P2 on the other hand has a finite non-trivial continuum 
limit = 0, Sq = 0. As for the coefficients t{ these are generally non-zero, 
since the diagram is a product of one-loop diagrams, but tQ = t\ = 0. For 
the leading coefficients we find 

rl = 0.000210(3), (4.4) 
s{ = 0.0003(14), (4.5) 
r{ = 0.003(9). (4.6) 

The extracted value of s{ above is unfortunately not very precise but consis- 
tent with the value 0.000646 obtained from the equation 

s{ + cf = (4.7) 

that we expect from 0(a) improvement. 

If we now analyze the series for + ci"^^*P2 with the known value for 
Ct^'^^* with the assumption that the equation ( [4.7| ) is fulfilled i.e. 

we get 

4^ = 0.000209(1), (4.9) 
= -0.0008(2). (4.10) 

Thus we obtain a refined estimate of P22 = Tq-' recorded in Table |^. Finally 
we estimate the improvement coefficient cp'^''* through 

^(2,2)* _ ifj ^ (^^ihi)*y ^ 0.0000(1) . (4.11) 
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4.1.2 The continuum behavior of p2i 

First of all, the expansion of pf has the same structure as that of i.e. it 
has a trivial continuum limit = 1 and also the coefficients = for all 
n. The analysis yields 

si = 0.0338(3), (4.12) 
rf = 0.115(2). (4.13) 

P2 on the other hand has a non-trivial finite continuum limit = 0, Sq = 
and = for all n. Our results for the leading terms are 

= -0.054641(1), (4.14) 
= -0.0252(2), (4.15) 
= 0.032(1). (4.16) 

The mass insertion term p^ has a linear divergence, n* = — 1 but s'_i = Sq = 
and fn = for all n. For the leading terms we find 

= 0.009568(1), (4.17) 
= 0.01198(11). (4.18) 

Finally the contribution p| has an expansion of the form 

pl = r^J + rl + sl ln(/) + ""i + 'i^^^) + . . . (4.19) 

t\ is expected to be zero because of tree level improvement and the data is 
consistent with this expectation. For this coefficient our analysis gives 

= 0.002586(2), (4.20) 
si = -0.0012(3), (4.21) 
rl = 0.0141(16). (4.22) 

On the other hand we know that the linear divergences must cancel and 
hence 

rt^ = -m^^Vii , (4.23) 



should hold. Indeed from ( 4.23 ) we find the value 0.0025841(3) which serves 
as an additional consistency check. In the next step we thus assume the 
cancellation of linear divergences and analyze 
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From this series we could get the still somewhat rough estimates of the con- 
tinuum behavior 

si' = -0.00104(12), (4.25) 
rg^ = 0.0099(7). (4.26) 

The Callan-Symanzik equation (|2.21|) (which incidentally requires to = 0) 
requires the logarithmic divergence s'^ to coincide with 26ii = —0.0010159, 
and our estimate ( [4.25| ) agrees with this within errors. We are thus justified 
to assume this to actually be the case and continue to analyze 

A + r^^A - HI) = rt + T + ^r-'W+'r-'^W + . . . (4,27) 

to obtain an improved estimate for the coefficient 

r^^ = 0.00978(2) . (4.28) 

The value for the coefficient p2i in Table |^ comes from 

P21 = r- + c«%^ (4.29) 

We also obtain (unfortunately rather poor) estimates for the 0(a) coefficients 

sf = -0.0016(17), (4.30) 
rf = -0.003(7). (4.31) 

In fact improvement requires 

- .r = ci'''^*4 + ci'''^'s{ + ci'>si + c«*.^ . (4.32) 



Using the previously obtained pf values = 0.16831(84) and s\ = 0.27848(40) 
gives sf" = —0.00103(7) which is consistent with but much more accurate 
than our estimate above. 

Finally the most accurate value we could obtain for improvement coeffi- 
cient c^^ was extracted by forming the combination 

P2 + m^c^P2 - 2&11 ln(/) + c['''^*p', + ci'''^*pl + d['^*pi + ci^*p', , (4.33) 

and analyzing the series with the assumption it has the form 

^ rf + snn(I) +tnn'^(I) 
ro + ^+ ' ' I — + ■■■ 4.34 

13 



thereby obtaining the estimate 

r^^ = 0.011(2) , (4.35) 

and finally 

^(2,1)* ^ 2c^°^*cf + \\ = 0.002(1) . (4.36) 
4.2 The case eg) = 

The analysis of the data for the case = is completely analogous to that 
for c^^ = 1 above, except for the fact that we can simply ignore also all other 
improvement coefficients. Thus e.g. for P22 only P2 has to be analyzed, and 
with the same ansatz as before we find 

= 0.000211(3). (4.37) 

For P21 we only take into account Ph- Again we find that the linear 

divergences cancel (with the appropriate value of m^^^) and the coefficient of 
ln(/) reproduces 26ii within errors. Finally we fit 

Pl + m^M - 26n ln(/) = r^, + + slHl)+ t^ln^I) ^ _ _ _ ^ ^^ ^^^ 

and the ensuing estimate of Tq gives us the value of P21 quoted in Table || for 
this case. 



5 Applications 

5.1 The relation of a-^ to asF 

To obtain the coefficients Cj(s) in eg. ( |1.1|) from our computation above, we 
have to use the known results for the coefficients appearing in the relation 
between and ao: 

ay^^{s/a) = ao + di{s)al + d2{s)al + . . . (5.1) 

with 

1 

diis) = ^iVf^'{-87r6orln(s) + rfi,}, (5.2) 

r=0 

1 

^2(5) = di{sf + Y.N/{-327i%irln{s) + d2r}. (5.3) 

r=0 
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The coefficients dir (and the /3— function coefficients) are given in Appendix 
A. 

The coefficients Ci{s) in eq.( |l.l| ) are now obtained through 

1 

ci(s) = -87rboHs)+Y.N/{dir-^mr}, (5.4) 

r=0 

2 

C2(S) = Ciisf -32TT\ln{s) + Y.N/{d2r-l67l^p2r}- (5.5) 

r=0 

Now it is clear that since eq. ( |1.1|) is a relation between continuum quantities 
all coefficients Cr must be independent of the lattice bare parameters, which 
serves as a further consistency check on the computations. This has already 
been observed at the 1-loop level i.e. from Table ^ and eqns.( |A^ , |A.8D one 
sees that du—Anpu is independent of At the two loop level, since ^22 = 
(there is no term oc A^f^ in ( |5.3|) ), we should find that ^22 is independent of Cgw 
Indeed our numerical results in Table ^ are consistent with this expectation. 
Finally ^21 — 167r^p2i should be independent of Csw; unfortunately this cannot 
be checked at present because the computation of the ^21 for Csw 7^ is not 
yet complete [^. However demanding the equality would require 



c^2iU=3,40)^, = 1.685(9) - 8.6286(2)c(;) , (5.6) 

which will serve as a good consistency check for the (general SU(A^)) compu- 
tation above in progress pSf. In fact the coefficient of above, 167r^rQ, is 



already independently checked because Tq should be related to the coefficient 
Ki in ( lOD through 

1 dKi(x) 



4:71 dx 



(5.7) 



x=l 

and this is fulfilled numerically to good precision. 

Putting all our numerical results together we find for A^ 



ci(s) = -87r6oln(s) + 1.255621(2) + 0.0398629(2)A'f, (5.8) 
C2(s) = ci(s)2 -3271^61 ln(s) 

+ 1.197(10) + 0.140(6)A'f - 0.0330(2)A'f2 . (5.9) 

This two-loop connection between the two different couplings determines the 
difference between the non-universal three-loop coefficients of their respective 
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/5-functions: 

Since the three-loop /3-function in the MS-scheme is known we can obtain the 
SF 3- loop beta function coefficient for the background field A (and 6 = ir/d) 
e.g. : 

&2kf=o = 0.482(7) X (47^)-^ (5.11) 
b2\Nt=2 = 0.064(10) X (47r)^3_ (5.12) 

The perturbative coefficients will find their application when the SF- 
coupling has been measured over a wide range of energies. Once the SF- 
coupling g is known for small box size, it is to be converted to the MS- 
coupling at high energy. One could conventionally choose the mass of the 
neutral weak boson Mz as a scale here. Another procedure — attractive for 
asymptotically free theories — is to extract the A-parameter which is simply 
related to the behavior at asymptotically large energy. It is a renormalization 
group invariant given by 



AsF = L-\bof)-''/'^''"^e-'/^'''~'^ exp |- J^' dg 



1 1 hi 
+ 



(5.13) 



The conversion to then amounts to an additional known factor: 



= 2.382035 (3) AsF, for A^f = 2. (5.15) 

If we insert a small g{L) belonging to a very small L (in physical units) into 
formula (|5.13|) , then the exponentiated integral is close to unity. Knowledge 



of the three loop term will give some estimate of the systematic error which 
will definitely be less than the statistical error in the near future. 

5.2 The step scaling function 

A central quantity in the ALPHA collaboration's approach is the step scaling 
function which generalizes the /9-function to finite rescalings, 

cr(s,n) = (?2(sL)|g2(i)=„. (5.16) 
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We remind the reader that we use a mass independent renormahzation scheme 
and set the quark mass to zero. On the lattice, a emerges as the continuum 
hmit of a finite lattice spacing approximant S, 

cr(s, m) = lim S(s, u, a/L). (5-17) 

Due to the absence of chiral symmetry before the continuum limit is taken one 
has to precisely specify a zero mass condition with the cutoff in place. In the 
numerical simulations a certain unambiguous definition based on the PCAC 
relation is adopted. This leads to / = L/a dependent expansion coefficients 
vrif^{l\ which only as / ^ oo go over to the convention independent values 



m, 



(0) 



and m^}^ in Table |T[ Only the latter are presently known to us. 
For the extrapolations in the previous section (including 0(a) improvement) 
this is of no concern. The following perturbative estimation of the full (all 
orders in a) lattice artifacts would however be more realistic with finite / 
mass expansion coefficients and is hence only given to get a first idea here. 
An improved version will be published elsewhere. 

All perturbative information about the convergence speed of S for s = 2 
is conveniently summarized in coefficients 5nj{a/L) defined by 

6{u,a L) = -— 5.18) 

n 

= E«"E^/^n,(a/L). (5.19) 

n=l j=0 

Setting Apij{I) = Pij{21) — Pij{I) we have 

5io = Apio-26ooln(2), 
6n = Apn-26oiln(2), 

520 = Ap2o-26ioln(2) -2Apio(Pio + &ooln(2)), 

521 = Ap2i -26nln(2) -2Apn(pio + 6ooln(2)) 

-2Apio(pii + 6oi ln(2)), 

522 = Ap22 -2Apn(pii + 6oiln(2)). 



\2 



(5 


20) 


(5 


21) 


(5 


22) 


(5 


23) 


(5 


24) 


in 


the 



All the Snj decay with an asymptotic rate proportional to (c _ 
0(a) improved theory. We are free to use S to cancel perturbative finite a 
effects from Monte Carlo data, as pointed out in P^|. For this purpose one 
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uses the same improvement coefficients as in the simulation and the expan- 
sion of rric for finite lattice spacing. Only for cases where nonperturbative 
improvement coefficients are used in the simulation, like Csw, one replaces 
these by the leading perturbative expressions. Fig.|I| gives an impression of 

L/a 

inf. 16 11 9 8 7 6 5 4 

% '■ 



X 

X 

X 

+ 

+ Nf=0 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 

{a/Lf 

Figure 1: Two loop contribution to lattice artifacts. 

the two loop artifacts for zero and two flavours. They are reasonably small 
and show the expected decay pattern. 
Acknowledgements 
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A. The 2-loop relation between aj^ and ao 

The coefficients dir appearing in the the relation between ajj^ and ao Eqs.( |5.2|) , 
(|5.3| ) have the form for general gauge group SU(A^) 

dn = K,, (A.2) 
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37r2 

d2o = ^ + k2 + hN\ (A.3) 

d2i = ^ + K,N. (A.4) 

The coefficients ki which have been computed in P| and checked by the 
authors of ref. |jlO| are given by 



ki = 2.135730074078457(2), (A.5) 
k2 = -2.8626215972(6), (A.6) 
^3 = 1.24911585(3). (A.7) 

The coefficient Ki is a function of c^^ and Wilson's parameter r. In our 
computations we always set r = 1 and hence only indicate the dependence 
on c(°) i.e. = Ki{4^) . We have ^ 

Ki{x) = -0.08414443(8) + 0.063419(l)x - 0.375024(1)^^ . (A.8) 

Similarly we have Kr = Kr{c'^^, c^^^) for r = 2, 3. So far only the pure Wilson 
case has been computed in ref. 

7^2(0,0) = 0.1890(2), (A.9) 
7^3(0,0) = -0.1579(3). (A.IO) 



The extension of the computation to the case Cs„ 7^ is now under way |28 
The universal coefficients of the beta-function are given by 

bi = ko + baNi, i = 0,l, (A.ll) 

with 

1 UN 



(A.12) 



(47r)2 3 ' 
^ 2 (A.13) 



(47r)2 3 

^To give a standard reference to the computation of Ki is difficult. The precise result 
in (|A. comes from Stefan Sint |2^. The result in eq.(29) of the recent paper |^ agrees 
with this . The computation of i^i(O) was first presented in refs. |^ and a more 

precise value is given in \ ^a\ . The extension of the calculation to Cgw 7^ was first presented 
in , however there are some errors in the result quoted there (of which the authors were 
conscious but unfortunately they did not publish an erratum). 
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1 S'lA'^ 

1 /13iV In , 

and the 3-loop beta function coefficient in the MS-scheme is given by 



with 



^Ms _ 1 2857iV3 

- (4vr)6 54 ' ^^-^^^ 
fWs _ 1 /1709iV^ 187 1 N 



^^22 - (4^)6 1 27 18iV^ ^^-^^^ 



B. The free fermion propagator in the abehan 
background field 

The Wilson-Dirac operator D, which is now taken to include the Sheikoleslami- 
Wohlert term, is diagonal in momentum and color space for the background 
field under consideration. We hence assume these fixed and discuss the in- 
version in time and Dirac-spin space. The calculational technique is adapted 
from the calculation the background field gluon propagator in |^ p!^ . 

The propagator will be entirely constructed from solutions to the homo- 
geneous equations 

(DV^/)(t) = with P+V^ (0) = (0) = 0, (B.l) 
{Dilj^){t) = with P^iIj\T) = ijtiT) = (B.2) 

with projectors P± = |(1 ± 7o)- In each case there are two independent 
solutions, which are constructed by forward (backward) recurrence in t for 
ipf [ipb) starting from the t = {t = T) boundary. We assemble the two 
solutions into two columns inipj , ipb- Using pseudo hermiticity, D^^ = 75D75, 
the defining equation 

DS{t, t') = 5t,e, < t, t' < T , (B.3) 
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and the fact that in a matrix sense S must also be the left inverse to D we 
derive 

Of, j ^■nt)V^P\t')h, t<t' 

^[t,t)- I ^;>(t)^t^/(t')t^5 t>t'. ^^-^^ 

Here ^ is a two by two matrix acting on the index that labels independent 
solutions. Inspection of the propagator equation at t' = t ± 1 yields in 
addition 

s{t, t) = i;iit)v^\tyj, + v^t(i) v^V^(t)W5 (B.5) 

and 

^i{t)vtp'_{ty-f5 = tpt{t)v^4ji{tyj5, (b.6) 

<(t)V^ViWW5 = V^i(t)V<(t)W5. (B.7) 
To determine V we exploit (p.3| ) at t = t', 

1 = -^i{t-i)v^t{ty^, + ^i{t-i)v^^i{t)h, (B.8) 



Now V can be isolated with the help of ( |B.6| ), (p.7| ) and 

^/(t - l)t75^/ (t) - ijf{t)h,iji{t - 1) = , (B.9) 

whose left hand side is shown to be t-independent by the homogeneous Dirac 
equation and vanishes for t = 1. The final result can be written as 

{V^y = ^^(l)^75<(0). (B.IO) 



C. Tables of expansion coefficients 

In this appendix we list perturbative finite lattice input data that went into 
our analysis. Numbers have been truncated such that roundoff errors affect 
the last digit only. 
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/ 


v\ 







3 





01406000332624868 


-0 


0182407224873773 





02748117066751495 


4 





01838971613803814 


-0 


0190839262114088 





02912977487693198 


5 





0213244240582013 


-0 


0182173759776234 





0272378718296281 


6 





0240118520229137 


-0 


0171194898351334 





0248632085828385 


7 





0266888600910630 


-0 


016021619058412 





0226936034816730 


8 





029366920333954 


-0 


014993754483268 





0208214421618160 


9 





032033898489622 


-0 


014058323610653 





0192191961300145 


10 





034684945753779 


-0 


013217310231851 





0178425484742938 


11 





037320642901569 


-0 


012464162363595 





0166512309905875 


12 





039943407218950 


-0 


011789492952359 





015612173517774 


13 





042555797334290 


-0 


011183638952885 





014698916407548 


14 





045159985052289 


-0 


010637686960534 





013890391073651 


15 





04775768717127 


-0 


010143794730699 





013169772817157 


16 





05035023257788 


-0 


009695210528629 





012523543034029 


17 





05293864891550 


-0 


00928617606722 





011940757139333 


18 





05552373638933 


-0 


00891179536156 





011412482658586 


19 





05810612401759 


-0 


00856790398265 





010931369933803 


20 





06068631107488 


-0 


00825095161082 





010491323683996 


21 





06326469747537 


-0 


00795790134914 





010087250557429 


11 





06584160619894 


-0 


00768614535208 





009714863852334 


23 





06841730002735 


-0 


00743343490502 





009370531360667 


24 





07099199418417 


-0 


00719782275317 





009051155905983 


25 





07356586599137 


-0 


00697761558627 





008754080827932 


26 





07613906232516 


-0 


00677133484944 





008477014637251 


27 





0787117066394 


-0 


00657768434489 





008217970511475 


28 





0812838974813 


-0 


00639552336330 





007975217365077 


29 





0838557242252 


-0 


00622384432090 





007747240012744 


30 





0864272578552 


-0 


00606175407446 





007532706527406 


31 





0889985593460 


-0 


00590845824714 





007330441330203 


32 





0915696803402 


-0 


00576324802706 





007139402877142 



Table 4: List of /-dependent coefficients -p^^ P2 and for c^^ — 1. 
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I 


P2 


P2 


P2 


3 


-0 


0472978469844831 


0.02748117066751495 


0.00036913190875759 


4 


-0 


0522497820522355 


0.0406644614491804 


-0.00019279262914480 


5 


-0 


0544132909735601 


0.0518711757366513 


0.000062922280309919 


6 


-0 


0555160735660078 


0.0625335624963019 


0.00024693874361661 


7 


-0 


0561252173383393 


0.0729651117622647 


0.00032829652675201 


8 


-0 


0564756935876461 


0.0832353491549959 


0.00035889080255099 


9 


-0 


0566798625246374 


0.0933761748042426 


0.00036802392907154 


10 


-0 


0567967356276605 


0.103414110425804 


0.00036829902719101 


11 


-0 


0568594539035442 


0.113371699857364 


0.00036493229776592 


12 


-0 


056887606070495 


0.123266779814805 


0.00036013282273909 


13 


-0 


056893249514880 


0.133112866383861 


0.00035486899918541 


14 


-0 


056884053136243 


0.142920056704081 


0.00034957573793481 


15 


-0 


056865026658050 


0.152695918495843 


0.00034444805110225 


16 


-0 


056839514713774 


0.162446191871116 


0.00033956882844280 


17 


-0 


056809789874931 


0.172175297258459 


0.0003349670866427 


18 


-0 


056777417868152 


0.181886691902227 


0.0003306456480143 


19 


-0 


056743488785825 


0.191583119155590 


0.0003265947619586 


20 


-0 


056708767026702 


0.201266784482244 


0.0003227989929269 


21 


-0 


056673790607995 


0.210939481638856 


0.0003192407628410 


22 


-0 


056638938178912 


0.220602684697517 


0.0003159021777408 


23 


-0 


056604474988832 


0.230257616284550 


0.0003127659514888 


24 


-0 


056570584880899 


0.239905298968026 


0.0003098158454325 


25 


-0 


056537392847533 


0.249546594496119 


0.0003070368462269 


26 


-0 


056504981113718 


0.259182234133082 


0.0003044152025451 


27 


-0 


056473400720539 


0.268812842375392 


0.0003019383874879 


28 


-0 


056442679941137 


0.278438955679912 


0.0002995950240854 


29 


-0 


056412830441785 


0.288061037388928 


0.0002973747949003 


30 


-0 


056383851821358 


0.297679489724326 


0.0002952683474247 


31 


-0 


056355734973936 


0.307294663501017 


0.0002932672016335 


32 


-0 


056328464590232 


0.31690686604973 


0.0002913636629602 



Table 5: List of /-dependent coefficients p\ and for — 1. 
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/ 


v\ 




P2 


3 





002930793635704460 





01113927564510310 


-0.00081995809255674 


4 





00648721649586021 





02176506549404950 


-0.00095088311404386 


5 





01042316816041390 





0326987791899857 


-0.00047704024765898 


6 





0147884221405374 





0441250945952050 


-0.00017002734632304 


7 





0194100790766652 





0557968441900837 


-0.000012571985852322 


8 





0241548411409164 





0675196648718815 


0.000069240759237721 


9 





028948291534668 





0791944332878618 


0.00011551757360719 


10 





033752733950337 





0907802022873036 


0.00014436115981723 


11 





038549666950648 





102264260790802 


0.00016392305375706 


12 





043330159467278 





113646262219098 


0.00017808308496443 


13 





048090083743491 





124930880299287 


0.00018883175180516 


14 





05282779555547 





136124585602541 


0.00019727388656343 


15 





05754297570840 





147234255392045 


0.00020407010116580 


16 





06223602393897 





158266580822338 


0.00020964219355829 


17 





06690772585583 





169227824652203 


0.00021427470438615 


18 





07155906252343 





180123737178240 


0.00021816831762828 


19 





07619109889499 





190959545294390 


0.00022146967768869 


20 





08080491792206 





201739975291261 


0.00022428900150278 


21 





08540158202310 





212469290287585 


0.00022671102155128 


22 





08998211127935 





223151332698423 


0.00022880209115049 


23 





09454747195837 





233789566873300 


0.00023061497375955 


24 





0990985714169 





244387119503159 


0.00023219217419543 


25 





1036362569094 





254946816719166 


0.0002335683174515 


26 





1081613167382 





265471217523167 


0.0002347718854825 


27 





1126744827534 





275962643571370 


0.0002358265093730 


28 





1171764335743 





286423205525826 


0.0002367519464219 


29 





1216677981308 





296854826274766 


0.0002375648294884 


30 





1261491592962 





307259261350171 


0.0002382792488587 


31 





1306210574330 





31763811686605 


0.0002389072090909 


32 





1350839937928 





32799286527967 


0.0002394589912651 



Table 6: List of /-dependent coefficients p^, pi and P2 for c^^ — 0. 
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D. Extrapolation of lattice perturbation theory 



In this appendix we discuss a modified method to extract the continuum 
hmit behavior of lattice Feynman diagrams. That is we assume to have 
precise results F{I) (for a given diagram or a sum of diagrams) for a range 
of I = L/a, Ii < I2 <...< In (typical subsets of I are 2, . . . , 32), and the 
goal is to determine the leading and sub leading behavior as J — > cxd. 

The precision of the data is limited by roundoff effects. We treat the 
roundoff errors Sf{I) as normally distributed superimposed noise, indepen- 
dent for different /. The most optimistic assumption for results from double 
precision (real*8) arithmetic would be 

5F{I) = e\FiI)\, e~10-l^ (D.l) 

but in general this is an underestimation since the evaluation of the diagrams 
involve large sums of terms of different signs. A comparison of double versus 
single (or quadruple) precision for a range of I can be used to derive a more 
realistic /— dependent e. (In our particular case we estimated a growth oc P 
for most contributions). 

For the extrapolation we assume an asymptotic expansion in functions 
|/fc(/) ,k = 1, 2, . . . , njj where successive terms become "smaller" with in- 
creasing index k, 

"/ 

F{I) = J2aMl) + R{I), (D.2) 

k=l 

where |i?(/)//„^. | ^ as / ^ 00. Of course, we have in mind something like 
{/i(^),/2(/),---} = l,ln///, ...,ln^J/J^, ...}. Denoting the n data values 
for F{I) as an n-dimensional column vector, the above equation reads 

F = fa + R, (D.3) 

where / (evaluated for the n values of / for which we have data) is regarded 
as an n X matrix and a is an n^-dimensional column vector that we want 
to determine. 



In ref. I^B] a recursive blocking technique was proposed to determine a 
which was claimed superior to making least square fits. In the following we 
shall find that blocking can actually be considered as a particular way of 
fitting with generalizations however being potentially more convenient. 
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We determine a by minimizing a quadratic form in the residues 



x' = {F-fa)'W\F-fa). (D.4) 

Here W'^ is an n-dimensional matrix of positive weights. It can be used to 
emphasize residues at small or large I in the minimization. Small / are less 
affected by roundoff, leading to less roundoff errors induced in a. On the 
other hand for larger / the asymptotic expansion holds to a better degree, 
and hence here smaller systematic errors are expected. In fact for our data 
we found little advantage in setting Wu = SjjP, —S<z<3; and thus we 
finally used a flat weight W = 1 and just changed the / range analyzed to 
observe convergence by moving /min while always keeping /max at the highest 
available size. The freedom in choice of W could however be useful in different 
situations and hence we keep it general in the following discussion. 
Minimization of (|D.4| ) leads to the equation 



fW'^fa = fW'^F. (D.5) 

We assume that the columns of Wf are linearly independent thus spanning an 
ra j-dimensional subspace [uf < n must always hold). Let P be the projector 
on to this subspace. Then 

Wfa = PWF, (D.6) 

is an equivalent equation that fixes a. A very stable and convenient way to 
solve for a is to construct the singular value decomposition for Wf . In 
a very simple version it amounts to a factorization 

Wf = USV^ , (D.7) 

where t/ is a column-orthonormal n x Uf matrix obeying 

U^U = 1; UU^ = P, (D.8) 

S is diagonal and V orthonormal, both of size nj x Uf. Then we get 

a = VS~^U^WF . (D.9) 

By simple error propagation the roundoff error of a is 
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A simple property of ( p.6| ) is the following. Imagine changing F by a 
component proportional to one of the functions fk included in the fit. Then 
it is easy to see that only the corresponding ak changes. In other words 
an estimate of ai from ( p.9|) is only uncertain due to components beyond 
/i, . . . , contained in F. Then it is clear that if we choose the /-window 
minimal {n = nj successive values) the solution is equivalent to a blocking 
procedure cancelling Uf — 1 components to isolate say ai. 

The estimation of the systematic errors is the most delicate problem in 
our context. After considerable experimentation we propose the following 
procedure. We assume that the remainder R can be modelled by a linear 
combination of the functions fuf+i, fnf+2, ■ ■ ■ , fuf+ur following those in- 
cluded in the fit. Typically in our case these would be the = 3 functions 
W ///™, In ///", 1//™, if 0(1//'"-^) is the last order in the fit. If they could 
be fitted individually (that is distinguished), they would be included in the 
analysis. Instead, for our error analysis we make separate fits including 
only one of them per fit in addition to the fi,i = 1, . . . ,nf. In this way we 
get the coefficients Ai, . . . , of the extra term in each case. The model re- 
mainders Axfrij^i-, A2fnf+2, ■ ■ ■ , ^n^/n^+n,. are fitted in the same way as the 
data ( p.9| ), and the maximal fit-coefficient of the Ur cases is taken as the 
systematic error da^^ of ak- 

We found that these errors based on the leading unincluded contribution 
are smaller in most (but not all) cases compared to the previous method 
focussing on the last included term. The present errors seem more realistic 
in the sense that expected relations (like coefficients of the /3-function) hold 
to an accuracy still conservatively but not grossly more accurate than the 
estimated errors. 
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